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We explore the phase diagram of neutral quark matter at high baryon density as a function of 
the temperature T and the strange quark mass M 3 . At T = 0, there is a sharp distinction between 
the insulating color-flavor locked (CFL) phase, which occurs where M 3 /fx < 2A, and the metallic 
gapless CFL phase, which occurs at larger / ll. Here, fx is the chemical potential for quark number 
and A is the gap in the CFL phase. We find this distinction blurred at T / 0, as the CFL phase 
undergoes an insulator-to-metal crossover when it is heated. We present an analytic treatment of 
this crossover. At higher temperatures, we map out the phase transition lines at which the gap 
parameters Ai, A2 and A3 describing ds-pairing, its-pairing and itd-pairing respectively, go to zero 
in an NJL model. For small values of Ms/ fx, we find that A2 vanishes first, then Ai, then A3. 
We find agreement with a previous Ginzburg-Landau analysis of the form of these transitions and 
find quantitative agreement with results obtained in full QCD at asymptotic density for ratios of 
coefficients in the Ginzburg-Landau potential. At larger Ms/fx, we find that Ai vanishes first, then 
A2, then A3. Hence, we find a "doubly critical" point in the (Mf /fx, T)-plane at which two lines of 
second order phase transitions (Ai — > and A2 — > 0) cross. Because we do not make any small-M s 
approximation, if we choose a relatively strong coupling leading to large gap parameters, we are able 
to pursue the analysis of the phase diagram all the way up to such large values of M„ that there are 
no strange quarks present. 

PACS numbers: 12.38.-t, 25.75.Nq 



I. INTRODUCTION 

Our understanding of the properties of matter at high 
baryon number density has developed rapidly in the last 
several years. At any densities that are high enough that 
nucleons are crushed into quark matter, the quark mat- 
ter that results must, at sufficiently low temperatures, 
be in one of a family of color superconducting phases . 
The essence of color superconductivity is quark pairing, 
driven by the BCS mechanism which operates whenever 
there is an attractive interaction between fermions at 
a Fermi surface. The QCD quark-quark interaction is 
strong and is attractive between quarks that are antisym- 
metric in color, so we expect cold dense quark matter to 
generically exhibit color superconductivity. Color super- 
conducting quark matter may well occur in the cores of 
compact stars. 

It is by now well-established that at asymptotic den- 
sities, where the up, down and strange quarks can be 
treated on an equal footing and the disruptive effects of 
the strange quark mass can be neglected, quark matter is 
in the color-flavor locked (CFL) phase, in which quarks of 
all three colors and all three flavors form Cooper pairs . 
The CFL phase is a color superconductor but is an elec- 
tromagnetic insulator, with zero electron density. 

To describe quark matter as may exist in the cores 
of compact stars, we need consider quark chemical po- 
tentials fx of order 500 MeV at most, meaning that the 
strange quark mass M s must be included: it is expected 
to be density dependent, lying between the current mass 
~ 100 MeV and the vacuum constituent quark mass ~ 
500 MeV. In bulk matter, as is relevant for compact stars 



where we are interested in kilometer-scale volumes, we 
must furthermore require electromagnetic and color neu- 
trality |3, |4| (possibly via mixing of oppositely charged 
phases) and allow for equilibration under the weak inter- 
actions. All these factors work to pull apart the Fermi 
momenta of the cross-species pairing that characterizes 
color-flavor locking. At the highest densities, we expect 
CFL pairing, but as the density decreases the combina- 
tion of nonzero M s and the constraints of neutrality put 
greater and greater stress on cross-species pairing, and 
reduce the excitation energies of those fermionic quasi- 
particles whose excitation would serve to ease the stress 
by breaking pairs [f|. 

If we imagine beginning in the CFL phase at asymp- 
totic density, reducing the density, and assume that CFL 
pairing is disrupted by the heaviness of the strange quark 
before color superconducting quark matter is superseded 
by the hadronic phase, the next phase down in density is 
the gapless CFL phase In this phase, quarks of all 

three colors and all three flavors still form Cooper pairs, 
but there are regions of momentum space in which certain 
quarks do not succeed in pairing, and these regions are 
bounded by momenta at which certain fermionic quasi- 
particles are gapless. The material is an electromagnetic 
conductor, with a nonzero electron density. 

Our goal in this paper is to map the phase diagram 
of neutral dense quark matter at nonzero temperature, 
answering the question of what phases and phase transi- 
tions result when CFL or gCFL quark matter is heated. 
In the case of CFL quark matter with strange quark mass 
M s = 0, this question has been answered in Ref. pj. In 
this most symmetric setting, there is a single phase tran- 
sition at a temperature T^ Fh below which there is CFL 



2 



pairing and above which there is no pairing. In mean 
field theory, 



T, 



CFL 



= 23— A , 

7T 



(1) 



where Ao is the CFL gap parameter at T — 0, estimated 
to be of order 10 to 100 MeV Q. The enhancement of 
T c /Ao by a factor of 2 1 / 3 over the standard BCS value 
(which our results confirm) originates in the fact that in 
the CFL phase with M s = there are eight fermionic 
quasiparticles with gap A and one with gap 2A 0. 
As at any phase transition at which a superconducting 
order parameter melts, gauge field fluctuations that are 
neglected in mean field theory can elevate T c and render 
the transition first order || . Because the gauge cou pling 
in QCD is strong, these effects are significant El E3- 
They have been evaluated to date only at M s = |3.ll0|. 
Once M s ^ 0, the mean field analysis alone becomes 
rather involved and we shall therefore leave the inclusion 
of fluctuations to future work. 

To see why M s ^ results in an intricate phase di- 
agram at nonzero temperature, we must introduce the 
pairingansatz that we shall use throughout this pa- 
per usiini: 



A 1 e a ' 31 e abl +A 2 e a ' 32 e ab2 +A 3 e a ^e ab3 (2) 



Here ip" is a quark of color a — (r,g,b) and flavor 
a = (u, d, s); the condensate is a Lorentz scalar, antisym- 
metric in Dirac indices, antisymmetric in color (the chan- 
nel with the strongest attraction between quarks), and 
consequently antisymmetric in flavor. The gap parame- 
ters Ai, A2 and A3 describe down-strange, up-strange 
and up-down Cooper pairs, respectively. At T = 0, 
there is an insulator-metal transition at M^/fi ~ 2Ai, 
at which the CFL insulator with A3 ~ A2 = Ai is re- 
placed by the gapless CFL metal, which has A 3 > A 2 > 
Ai > [1 EL (If the CFL phase is augmented by a K°- 
condensate |l2( , the CFL— >gCFL transition is delayed to 
Ml) '// ~ 8Ai/3 0.) An analogous zero temperature 
metal insulator transition has been analyzed in Ref. . 

At M s = 0, Ai = A 2 = A 3 = A CFL , and A CF l 
decreases from Ao at T = to at T C CFL . As soon 
as M s 7^ 0, however, we can expect that Ai, A 2 and 
A3 do not all vanish at the same temperature. This 
expectation is evident for the gCFL phase, which has 
A 3 > A 2 > Ai > already at T = 0. We shall see that 
it also applies to the CFL phase with M s 7^ 0. 

The distinction between an insulator and a metal is 
sharp only at T = 0. At any nonzero temperature, 
there will be some nonzero density of thermally excited 
fermionic quasiparticles, some of which are charged. This 
means that the CFL^gCFL "transition" should be a 
crossover at any nonzero temperature 0]. The CFL 
phase at T = with M s ^ has fermionic quasiparticles 
with opposite charges whose excitation energies differ. 
This means that upon heating this phase, the chemical 
potentials needed to maintain neutrality are not the same 



as at zero temperature. Adjusting the chemical poten- 
tials feeds back into the gap equations for Ai, A 2 and A 3 
differently, and these gap parameters can therefore have 
different T-dependence. If we start in the CFL phase at 
T = 0, heat, then increase M^/fj, above that for the zero 
temperature CFL— >gCFL transition, and then cool back 
to T = we can go from CFL to gCFL without ever 
crossing a phase boundary. We illustrate this by showing 
an example phase diagram in Fig. ^ 

The purpose of this paper is to derive and understand 
the phase diagram of Fig. ^ We shall present the model 
and approximations that we use in Section [H] We fol- 
low the conventions for naming phases used in previous 
literature: 

A 1 ,A 2 ,A 37 ^0 (g)CFL 
Ai=0, A 2 ,A 37 ^0 uSC 
A 2 = 0, A 1; A 3 ^0 dSC 
Ai = A 2 = 0, A 3 + (g)2SC 
Ai = A 2 = A 3 = UQM, 

with "UQM" meaning unpaired quark matter. The ori- 
gin of the remaining names is that in the 2SC phase, only 
quarks of two flavors and two colors pair, whereas in the 
uSC (or dSC) phase, all u (or all d) quarks pair 
The phase diagram features three lines denoting second 
order phase transitions at which Ai, A 2 or A3 vanish. At 
M s = 0, all three gaps vanish at the same temperature 
but at generic values of M 3 , there are three distinct tran- 
sition temperatures. Two of these lines cross at a point 
denoted by the circle: at this "doubly critical point", Ai 
and A 2 vanish at the same temperature. The phase di- 
agram is intricate, and one natural question is to what 
extent its features are generic. We address this by vary- 
ing parameters, as we shall discuss in Sections llIII and lTvl 
We shall find that the physics at large Mf / fj,, where we 
find a first order phase transition in Fig.Q] is not generic, 
changing qualitatively for larger values of Ao- We find 
that the shape of the second order phase boundaries sep- 
arating the CFL, dSC, uSC and 2SC phases is generic, 
as is the existence of the doubly critical point. 

Physics in the vicinity of any of these phase transi- 
tions can be analyzed using a Ginzburg-Landau approx- 
imation, in which the relevant A or A's are taken to be 
small. This analysis has been performed at M s = in 
Refs. HHE3 and at small but nonzero M s in Ref. [l6| . 
We shall show in Section IVT1 that our numerical results 
for the three phase transition lines at small M 2 / fj,Ao are 
in quantitative agreement with the Ginzburg-Landau ap- 
proximation, and shall show that the ratios of coefficients 
in the Ginzburg-Landau potential that we obtain in our 
model agree quite well with those obtained in full QCD 
at higher densities. At small M 2 /fiAo, the region of the 
phase diagram where the Ginzburg-Landau analysis ap- 
plies near T c , A 2 vanishes at a lower critical temperature 
than Ai. In heating gCFL quark matter, however, we ex- 
pect Ai to vanish first because it is already much smaller 
at T = 0. Our model analysis can be extended to larger 
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FIG. 1: Phase diagram of dense neutral quark matter in the (Ms / [J,, T) plane. Only the solid curves represent phase transitions. 
The thin solid curves denote mean field second order phase transitions at which Ai or A2 or A3 vanish. Two of these lines 
cross at the circle, which we call a "doubly critical" point. The heavy solid line is a first order phase transition that ends at 
a tricritical point denoted by the diamond. The dashed curves indicate locations where, as Mg/fj, is increased, new gapless 
modes appear. At any nonzero temperature, however, there is no physical distinction between a phase with truly gapless modes 
and one in which the same modes have excitation energies of order T. The dashed curves therefore have physical significance 
only where they cross T = 0. In this phase diagram, this occurs only at the triangle, which denotes the zero temperature 
CFL^gCFL transition, which is a quantum critical point. The dotted line separates the Ai > A2 and A2 > Ai regions. This 
phase diagram is drawn for fi — 500 MeV, with M 3 and T varying, and with the quark-quark interaction strength fixed and 
chosen such that the CFL gap is Ao = 25 MeV at T = 0. We shall describe the model and approximations within which this 
phase diagram has been obtained in Section [H] and describe the results summarized by this diagram in Sections ITTT1 and HYI 



M s and smaller T than the Ginzburg-Landau analysis, 
allowing us to see how the phase diagram fits together, 
consistent with this expectation. 

A phase diagram similar to ours was obtained previ- 
ously in Ref. [15j , but our results differ in a crucial, qual- 
itative respect at low M^j\i, where we believe that our 
results are robust: the authors of Ref. 0] found that Ai 
always vanishes at a lower temperature than A2 , meaning 
that there is no dSC region on their phase diagram. This 
also disagrees with the Ginzburg-Landau analysis |47| . 
We shall detail our approximations in Section |TT1 but we 
note already here that whereas in Refs. P.l6lll5j| M„ was 
assumed to be much smaller than /x, we do not make 
such an approximation. This likely explains the differ- 



ences between our results and those of Ref. [15| at large 

Although not associated with a phase transition, it is 
interesting to ask how the CFL phase ceases to be an 
insulator as it is heated. This must happen, give n the 
phase diagram that we and the authors of Ref. [l5j find, 
but how? At small T, the number density of charged 
quark quasiparticles in the CFL phase is exponentially 
small. Because of the interplay between M s ^ and the 
constraints of neutrality, the excitation energies of oppo- 
sitely charged quasiparticles are not the same and these 
thermally excited quasiparticles have a net charge which 
must be balanced by a nonzero electron density of order 
yitg, where jjL e is the electron chemical potential. Because 
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the quasiparticle densities are proportional to the quark 
Fermi surface area ~ /i 2 , where > /j e , we show that fj, e 
ceases to be exponentially suppressed - rising rapidly to 
fi e ~ M 2 /4/i - in an insulator- metal crossover that occurs 
in a narrow, and quite low, range of temperatures. We 
describe this crossover analytically in Section There 
are other charged excitations in the CFL phase, namely 
the pseudo-Nambu-Goldstone bosons, that are excited 
when the system is heated. We argue quantitatively that 
their effects are small. 

In Section[n]we detail our model and approximations, 
stressing also those approximations used in previous work 
that we have not made. In Section II I II we review re- 
sults at T = 0. In Section IIVI we present our results 
at T 7^ 0, analyzing the phase diagram in Fig. ^ an d 
that for other values of A . In Section we analyze 
the insulator-metal crossover that occurs when the CFL 
phase is heated. In Section IVII we make the connec- 
tion between our work and the Ginzburg-Landau analysis 
quantitative. We conclude in Section IVTT1 



II. MODEL AND APPROXIMATIONS 

A. Pairing ansatz for gap parameters 

We assume that the predominant diquark condensate 
is a Lorentz scalar (antisymmetric in Dirac indices), is 
antisymmetric in color (as is favored energetically), and 
consequently is antisymmetric in flavor. That is, we as- 
sume a condensate of the form PJl. The gap parameters 
Ai, A2 and A3 describe a 9 x 9 matrix in color-flavor 
space that, in the basis (ru, gd,bs,rd, gu,rs,bu, gs,bd), 
takes the form 
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(3) 

We see that (gs,bd), (rs,bu) and (rd,gu) quarks pair 
with gap parameters Ai, A 2 and A3 respectively, while 
the (ru, gd, bs) quarks pair among each other involving 
all the A's. The most important physics that we are 
leaving out by making this ansatz is pairing in which 
the Cooper pairs are symmetric in color, and therefore 
also in flavor. Diquark condensates of this form break 
no new symmetries, and therefore must arise in the CFL 
phase However because the QCD interaction is re- 

pulsive between quarks that are symmetric in color these 
condensates are numerically insignificant 0, ITEl llgj . 



B. CFL symmetries and excitations 

To set the stage for our analysis, and to introduce the 
excitations that shall concern us, we briefly summarize 
the properties of the CFL phase 0. If we set all three 
quark masses to zero, the diquark condensate in the CFL 
phase (which then has Ai = A 2 = A3) spontaneously 
breaks the full symmetry group of QCD, 

[S*7(3)coior] x SU(3) L x SU(3)h x U(1) b 

^™ (4 ) 

3 P(1)q} 

where SU(3) co \ OT and electromagnetism U(1)q are 
gauged, and the unbroken SU(3) c +l+r subgroup con- 
sists of flavor rotations of the left and right quarks with 
equal and opposite color rotations, and contains an un- 
broken gauged "rotated electromagnetism" U(1)q [3.ll9|. 
The CFL phase has the largest possible unbroken sym- 
metry consistent with diquark condensation, achieved by 
having all nine quarks participate equally in the pairing, 
and this gives the maximal pairing free energy benefit. 
Not surprisingly, ab initio calculations valid at asymp- 
totic densities confirm that the CFL phase is the ground 
state of QCD in the high density limit fTL l2p| . 

The axial flavor U(1)a rotations are not a symmetry 
of QCD, as they are explicitly broken at all densities by 
instanton effects. At asymptotic densities, these effects 
become small but at accessible densities they cannot be 
neglected H El El and they shall play an im- 

portant role in Section IV AI 

In the limit of three massless quarks described above 
there are 17 broken symmetry generators in the CFL 
phase, 8 of which become longitudinal components of 
massive gauge bosons and 9 of which remain as Nambu- 
Goldstone bosons. The massive gauge bosons have ex- 
citation energies of order gfi, where the QCD gauge 
coupling g is not small, and we neglect them through- 
out. In the real world, with its two light quark fla- 
vors with masses < 10 MeV, and a medium-weight fla- 
vor, the strange quark, with mass > 100 MeV, the 
eight Nambu-Goldstone bosons coming from the break- 
ing of chiral symmetries acquire masses 0, 0, 0, 
l25l |2(| , and furthermore the CFL condensate may ro- 
tate in the /^-direction within the manifold describing 
these mesons [IJ- We neglect the possibility of K°- 
condensation throughout. We shall see in Section IV AI 
that this corresponds to assuming that the instanton con- 
tribution to the K° mass is comparable to Ao E3 ! 
which is likely the case for all but the largest Ao that 
we consider. We must also consider the thermally ex- 
cited pseudo-Nambu-Goldstone bosons, but we show in 
Section IV AI that their effects are negligible. The ninth 
Nambu-Goldstone boson, corresponding to the sponta- 
neous breaking of U(1)b and hence to superfluidity, re- 
mains massless even once quark masses are taken into 
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account and therefore plays a crucial role in many low 
energy properties of the CFL phase, for example in its 
specific heat, thermal conductivity, neutrino opacity, and 
neutrino emissivity at low temperatures [2?1 I2M |2t|. 
However, because it is neutral it will play no role in our 
analysis. 

The excitations of the CFL phase that are most rele- 
vant to our considerations are the nine fermionic quasi- 
particles, whose excitation requires the breaking of pairs. 
In analyzing the response of the CFL phase to the strange 
quark mass at zero temperature, the quasiquark disper- 
sion relations signal the instability corresponding to the 
disruption ofpairing that constitutes the CFL— >gCFL 
transition 5,6]. These dispersion relations are even more 
crucial at nonzero temperature, because they determine 
the number densities of thermally excited quasiparticlcs, 
some of which are charged and which we shall see there- 
fore play a crucial role in the neutrality conditions. Fur- 
thermore, these dispersion relations define the gap equa- 
tions that we shall solve at nonzero temperature, again 
as we shall see. 



C. Color and electric neutrality 

Stable bulk matter must be neutral under all gauged 
charges, whether they are spontaneously broken or not. 
Otherwise, the net charge density would create large elec- 
tric fields, making the energy non-extensive. In the case 
of the electromagnetic gauge symmetry, this simply re- 
quires zero charge density. In the case of the color gauge 
symmetry, the formal requirement is that a chunk of 
quark matter should be a color singlet, i.e., its wave- 
function should be invariant under a general color gauge 
transformation. Color neutrality, meaning equality in the 
numbers of red, green and blue quarks, is a less stringent 
constraint. A color singlet state is also color neutral, 
whereas the opposite is not necessarily true. However 
it has been shown that the projection of a color neutral 
state onto a color singlet state costs no extra free energy 
in the thermodynamic limit 30]. Analyzing the conse- 
quences of the requirement of color neutrality therefore 
suffices for our purposes. 

In nature, electric and color neutrality are enforced 
by the dynamics of the electromagnetic and QCD gauge 
fields, whose zeroth components serve as chemical poten- 
tials which take on values that enforce neutrality 0, . 
Since we are limiting ourselves to color neutrality and not 
color singletness we have to consider only the U(l) x U(l) 
diagonal subgroup of the color gauge group. This sub- 
group is generated by the diagonal generators T 3 = 
diag(i,-i,0) and T 8 = diag(±,±,-§) of the SU(3) 
gauge group. Electromagnetism is generated by Q = 
diag(|, — -|, — |) in flavor space (u, d, s). The zeroth com- 
ponents of the respective gauge fields serve as chemical 
potentials fj,^ and /x 8 coupled to T3 and T 8 charges, and as 
an electrostatic potential /i e coupled to the negative elec- 
tric charge Q. (We make this last choice so that jx e > 



corresponds to a density of electrons, not positrons.) The 
dynamics of the gauge potentials then require that the 
charge densities, which are the derivatives of the free en- 
ergy with respect to the chemical potentials, must vanish: 



n 8 



(5) 



In an NJL model with fermions but no gauge fields that 
we shall employ, one has to introduce the chemical po- 
tentials /i e , /i3 and /x 8 "by hand" in order to enforce color 
and electric neutrality in the same way that gauge field 
dynamics does in QCD Q. 

A generic diquark condensate will be neither electri- 
cally nor color neutral, so it will spontaneously break 
these gauge symmetries. However it may be neutral un- 
der a linear combination of electromagnetism and color. 
Indeed, any condensate of the form is neutral with 
respect to the "rotated electromagnetism" generated by 
Q = Q — T3 — TjTg, so U (1)q is never broken. This means 

that the corresponding gauge boson (the "Q-photon" ) , a 
mixture of the ordinary photon and one of the gluons, 
remains massless. In both the CFL and gCFL phases, 
the rest of the SU{?>) co \ OT x U(1)q gauge group is spon- 
taneously broken, meaning that the combination of the 
photon and gluons orthogonal to the Q-photon, and all 
the other gluons, become massive by the Higgs mecha- 
nism. 

In the CFL phase at T — 0, many consequences of 
the neutrality equations can be analyzed in a model- 
independent way A central result is that the free 
energy Q is independent of 



\{jl e +/X 3 + § /is) I 



(6) 



the chemical potential for Q-charge, as long as fig is not 
so large that it disrupts CFL pairing. This means that 
the CFL phase is a Q-insulator at T = [3J| . Translating 
back into fi e , ^3 and /x 8 , this means that two of these 
three can be fixed as functions of the third: 



M3 = Me, 



Ms = 



Ml 

2jJL 



Me 
2 



(7) 
(8) 



where ^ e can lie anywhere within the range — A2 + -grj- < 
He < +A3. If the lower (upper) bound of this range is 
violated rs-bu (rd-gu) pairing is disrupted: this range can 
be thought of as the bandgap for the CFL insulator 
The free energy D, has a plateau within this range of /Iq . 
Electrons add a new term ~ /z^ to fl, adding a gentle 
curvature to the plateau, and selecting /j, e = 0, meaning 
zero electron density, as the correct solution for neutral 
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CFL matter at T = 0. (/X3 and /is are then given by 
(|7I8|) with fj, e = 0.) At finite temperature where quasi- 
quarks with non-trivial Q are excited, we shall see that 
because fl has an only- very-gently-curved plateau, it will 
not take much thermal excitation of quasiquarks to shift 
fi e significantly away from zero. 



D. A Model for the Thermodynamic Potential 

We are interested in physics at non-asymptotic den- 
sities, and therefore cannot use weak-coupling methods. 
We are interested in physics at zero temperature and high 
density, at which the fermion sign problem is acute and 
the current methods of lattice QCD can therefore not 
be employed. For this reason, we need to introduce a 
model in which the interaction between quarks is simpli- 
fied, while still respecting the symmetries of QCD, and 
in which the effects of M s , /x e , fi^, /xg and T on CFL 
pairing can all be investigated. The natural choice is to 
model the interactions between quarks using a point-like 
four-fermion interaction, which we shall take to have the 
quantum numbers of single-gluon exchange. We work 
in Euclidean space. Our partition function Z and free 
energy density £1 are then defined by 

Z = c- f3vn =Af J V^jVifiexp^J C(x)d 4 



(9) 



c(x) = ^(i0 + |i-M)v-gG$r^)$r^) 



where the fields live in a box of volume V and Euclidean 
time length (3 — l/T, and i/l = /i.74. The interaction 
vertex has the color, flavor, and spin structure of the 
QCD gluon-quark coupling, = r ) ll T A with the Ta 
normalized such that Tr TaTb = %&ab- The mass ma- 
trix M = diag(0,0, M s ) in flavor space. The chemical 
potential fi is a diagonal color-flavor matrix depending 
on /i, fi e , /13 and /is. The normalization of the four- 
fermion coupling 3G/8 is as in the first paper in Ref. 0. 
In real QCD, the ultraviolet modes decouple because of 
asymptotic freedom but, in the NJL model, we have to 
add this feature by hand, through a UV momentum cut- 
off A in the momentum integrals. The model therefore 
has two parameters, the four-fermion coupling G and the 
three-momentum cutoff A, but it is more useful to pa- 
rameterize the interaction in terms of a physical quan- 
tity, namely the CFL gap parameter at M a = T = at 
a reference chemical potential that we shall take to be 
500 MeV. We shall call this reference gap Ao. We have 
checked that if we vary the cutoff A by 20% while simul- 
taneously varying the bare coupling G so as to keep Ao 
fixed, then our results change by a few percent at most. 
All the results that we present are for A = 800 MeV. 
Most of the results that we present are for a coupling 
strength chosen such that Ao = 25 MeV, but we shall 
also discuss results obtained with stronger couplings for 
which A = 40 MeV and A = 100 MeV. 



The free energy ft obtained from the Lagrangian © 
upon making the ansatz @ for the diquark condensate 
and working in the mean field approximation has been 
presented in Ref. [I . More sophisticated derivations exist 
in the literature [l|, but since we are assuming that the 
only condensate is of the form J2J) we can Fierz transform 
the interaction to yield products of terms that appear in 
(J2J) , and discard all the other terms that arise in the Fierz 
transformed interaction which would anyway vanish after 
making the mean field approximation. This yields 



G 



(10) 



where 



(Pn)?f = Cj 5 £ aPv e t]V (no sum over V ) (11) 

and P v = 74P^74. The index r\ labels the pairing channel: 
i] = 1, 2, and 3 correspond to d-s pairing, u-s pairing, and 
u-d pairing. The gap parameters can then be defined 
precisely as 

A n = ±G(V> T P^) . (12) 
As derived in Ref. , the free energy can then be written 



as 



n = --^ /Vp 2 E{l^(ril+ 2Tlll ( 1 + eH£3(?)l/r )} 

T JO j 
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G ' 



-(Al+Al+Al) 



12tt 2 



G 



180 



(13) 



where the electron contribution for T ^ has been in- 
cluded. The functions £j(p) are the dispersion relations 
for the fermionic quasiparticles. They are not explicitly 
T-dependent, but they do depend on the A's and //'s 
which are T-dependent. The quasiquark dispersion re- 
lations Ej(p) are the values of the energy at which the 
propagator diverges, 

dct S- 1 {-is jip)^) = . (14) 

Here, the inverse propagator S 1-1 is the 72 x 72 matrix 



P A 



(15) 



that acts on Nambu-Gor'kov spinors 



* = 



jp(p) 
4> T (~P) 



f = 



(tf(p) i> T (-p)) . (16) 



The sum in (|13|l is understood to run over 36 roots, be- 
cause the Nambu-Gor'kov formalism has artificially made 
each \Ej\ doubly degenerate. There is a further (physi- 
cal) degeneracy coming from spin, meaning that there 
are only 18 distinct dispersion relations to be found. The 
free energy of any solution to the gap equations obtained 
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from this mean field free energy is the same as that ob- 
tained from the Cornwall- Jackiw-Tomboulis effective po- 
tential j^|, with the mean field approximation imple- 
mented via dropping the contribution of 2PI diagrams. 

We shall evaluate f2 numerically, meaning that at every 
value of the momentum p we must find the zeros of the 
inverse pro pag ator. In order to do this numerically, we 
follow Ref. and note that it is equivalent and faster to 
find the eigenvalues of the "Dirac Hamiltonian density" 
H(p), which is related to S^ 1 by 



det S~ 1 (p i ,p) = det 7 4 (ip 4 - H{p)) 



(17) 



which makes it clear that eigenvalues of Ti(p) are zeros of 
S -1 . Like S^ 1 , the Dirac Hamiltonian H(p) is a 72 x 72 
matrix in color, flavor, spin, and Nambu-Gor'kov space. 
This matrix can be decomposed into 4 blocks - three 
16 x 16 matrices for the (rd-gu), (rs-bu) and (gs-bd) pair- 
ing sectors, and one 24 x 24 matrix for the (ru-gd-bs) sec- 
tor. The absolute values of the eigenvalues are quadruply 
degenerate due to the Nambu-Gor'kov doubling and the 
spin degeneracy. 

We have found considerably simpler expressions for the 
(rd-gu), (rs-bu) and (gs-bd) sectors. The calculation be- 
comes easier if one adopts a different choice of Nambu- 
Gor'kov basis, i.e. * = (ip(p), Cip T (-p)) T as in Ref. [34| . 
Then the Nambu-Gor'kov and spin degeneracies are man- 
ifest, as the 16 x 16 matrix separates into 4x4 blocks 
with the form 



H 



(gs-bd) 



■fJ-bd 


P 





iAi 


P 


— ^bd 


-iAi 








zAi 


flgS + M S 


P 


iAi 





p 


flgS- M S 



(18) 



for the (gs-bd) sector and likewise for other sectors. Each 
of the 4 eigenvalues of this matrix, and of its counterparts 
for the other sectors, contributes twice (for spin) in the 
sum in (|13fl . The quark chemical potentials occurring in 
(|T%)l and below are defined straightforwardly: 



jl ru 


= 


§Me + 


5^3 4 


■ 3^8, 


fJ-gd 


= M + 


3^e - 


5M3 4 


- 


Hbs 


= H + 


\fie - 


3^8, 




Mrd 


= M + 


§Me + 


3^3 4 




l^gu 


= M~ 




5^3 4 




l^rs 


= M + 


§Me + 


5^3 4 


■ 3M8, 


Hbu 


= M~ 




2 

3^8, 






= M + 


jMe - 


5^3 4 






= M + 


3^e — 


§M8 ■ 





(19) 



In the (ru-gd-bs) sector, we could not find any sim- 
ple way to make the spin degeneracy manifest, but the 
Nambu-Gor'kov degeneracy is manifest as the 24 x 24 
matrix separates into 12 x 12 blocks with the form 



\x ru 


P 





-iAa 





-iA 2 


p 


f^ru 


iA 3 





iA 2 








-iA 3 


[i-gd 


P 








iA 3 





P 


Hgd 











-iA a 








flbs+Ms 


p 


iA 2 











P 


fi bs - M s 











-iA l 














zAi 


























-iAi 














iAi 























Q 




































\ 

























-iAi 
















iAi 




























-iAi 
















iAi 













fl bs + M s 


p 











-iA 2 




P 


-fibs -M s 








iA 2 













-Pgd 


p 





-iA 3 










P 




iA 3 










-iA 2 





-iA 3 


f^ru 


P 




iA 2 





iA 3 





p 


\X TU 


) 



The 12 eigenvalues of this matrix each contribute once in 
the sum in (|13|l . They occur in degenerate pairs due to 
spin. 



A stable, neutral phase must minimize the free energy 
(|13|l with respect to variation of the three gap parameters 



Ai, A 2 , A3, meaning it must satisfy 

on on on 



9Ai 



0, 



dA 2 



0, 



dA a 







(21) 



and it must satisfy the three neutrality conditions 
The gap equations (|21|l and neutrality equations (jSJ) form 
a system of six coupled integral equations with unknowns 
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the three gap parameters and ^t e ,/i3 and /is- We have 
solved these equations numerically at a grid of values of 
T and M^/fi. We evaluate using <fT3*)>. finding the Ej 
by determining the eigenvalues of the matrices specified 
explicitly above. We evaluate partial derivatives of f2 us- 
ing finite difference methods. (We used three and five 
point finite difference evaluations of derivatives, with the 
interval in the A or fi with respect to which fl is being 
differentiated taken as 0.04 MeV.) We have worked en- 
tirely with /i = 500 MeV, varying M^/fi by varying M s , 
and have worked at three values of the coupling G chosen 
so that Ao, the gap parameter at M s = T = 0, takes on 
the values A = 25, 40 and 100 MeV. 

E. Approximations made and not made 

Having described our calculation in explicit detail, we 
close this section by enumerating the approximations 
that we are making, and stressing several that we have 
not made. 

We have made the following approximations: 

1. We work throughout in a mean field approximation. 
Gauge field fluctuations will surely be important, 
and must be included in future work. 

2. We neglect contributions to the condensates that 
are symmetric in color and flavor: these are known 
to be present and small, smaller than the contribu- 
tions that we include by at least an order of mag- 
nitude nmm. 

3. We choose light quark masses M u = Ma = and we 
treat the constituent strange quark mass M s as a 
parameter, rather than solving for an (ss) conden- 
sate. These approximations should be improved 
upon, along the lines of Refs. [ill l34j . In future 
studies in which the (qq) condensates are solved 
for dynamically, the six-quark 't Hooft interaction 
induced by instantons must be included, since it in- 
troduces terms of the form (ss) | A3 | 2 , for example. 

4. We ignore meson condensation in both the CFL 
phase and gCFL phases. We shall see in Sec- 
tion lV Al that in the CFL phase this means that we 
are assuming that the instanton contribution to the 
K° mass is comparable to or larger than Ao • 
Meson condensation in the gCFL phase has yet to 
be analyzed. 

We expect that these approximations have quantitative 
effects, but none preclude a qualitative understanding of 
the phase diagram. 

The authors of Refs. 0, |(| made further simplifying 
assumptions. They worked only to leading nontrivial 
order in Ai, A2, A3, /j, e , ^3 and /is- They neglected 
the effects of antiparticles. And, most serious, they in- 
corporated M s only via its leading effect, namely as an 
effective shift — M 2 //i in the chemical potentials of the 



strange quarks. This limits the regime of applicability 
of their results to M s < \i. For A ~ 100 MeV, there 
is pairing even at M s > fj, and our calculation runs into 
no difficulties in this regime. (We shall show some re- 
sults, even though this likely corresponds to such small fi 
that in reality the hadronic phase has taken over.) In the 
next section, we shall present our T = results. Those 
with A = 25 MeV, as in Refs. [H, , are in very good 
agreement with the results of Refs. This provides 

quantitative evidence for the validity of the approxima- 
tions made in Refs. H @ at A = 25 MeV. 



III. ZERO TEMPERATURE RESULTS 

In this section we shall present the results of our model 
analysis at T = 0. Here and throughout this paper, 
we set fj, = 500 MeV corresponding to a baryon density 
about ten times that in nuclear matter, which is within 
the range of expectations for the density at the center of 
compact stars. In reality, both M s and Ao vary with /x, 
and our knowledge of both is uncertain. Our model is by 
construction well-suited to study the disruptive effects of 
M s on pairing but it is not adequate to fix the density- 
dependent values of either M s or An quantitatively. In 
this study, therefore, as in Refs. |I|5j,|SlL5| we treat both 
M s and Ao simply as parameters. We shall plot our re- 
sults versus M 2 //i, rather than versus M s itself, because 
the strength of the disruptive effects introduced by the 
strange quark mass are characterized by M^/fi. For ex- 
ample, the splitting between Fermi surfaces in unpaired 
quark matter is of order M 2 /4/i, and the CFL^gCFL 
transition occurs where M 2 //z ~ 2Ai. Although we only 
present results at a single value of //, because the effects 
of M s are controlled by M 2 / fi, and because in reality M s 
is itself a decreasing function of /i, one can think of our 
plots as giving a qualitative description of the effects of 
varying /i, with high density at small M^/fi and low den- 
sity at large M 2 ///. This description is only qualitative 
because our plots are each made with some fixed value 
of Ao, whereas in reality Ao depends on fi. 

We begin by choosing the four fermion coupling G 
so that A = 25 MeV, as in Refs. d @. The solid 
curves in Fig. [3 show the gap parameters as a function 
of M s 2 //i, and those in Fig. [31 show the chemical poten- 
tials. Theseplots are in very good agreement with those 
of Refs. [Hill, indicating that the approximations made 
in that paper that we have dispensed with here, chiefly 
the small- Af 2 / assumption, were good approximations 
for A = 25 MeV. For small M 2 //x, we see the CFL 
phase with Ai = A2 ~ A3. The small difference be- 
tween Ai = A2 and A3, less than 2% everywhere within 
the CFL phase, is an example of an effect that we can 
see but that cannot be seen in the small- M 2 / /x 2 approx- 
imation in which M s is approximated as a shift in the 
chemical potential for the strange quarks @, The 
CFL^gCFL transition occurs at M 2 //i = 46.8 MeV, 
which is very close to 2Ai, since Ai = 23.2 MeV at 
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FIG. 2: Gap parameters Ai , A2 and A3 as a function of Ml / fx 
at T — 0. The solid curves show the CFL/gCFL solution, with 
the CFL— ^gCFL transition occurring where M%//j. ~ 2Ai. 
There is a first order phase transition between the gCFL phase 
and unpaired quark matter at Ml /fx = 125.3 MeV, denoted 
by the the thin vertical line. To the right of this line, the 
gCFL solution is metastable. We also find a neutral 2SC so- 
lution, with A3 given by the dashed curve in the figure, which 
undergoes a transition to the gapless 2SC phase of Refs. |3fi| 
at Ml In = 112 MeV. However, from Fig. we see that the 
(g)2SC solution is everywhere metastable, having a larger free 
energy than the (g)CFL solution at the same Ml / fx. 

Mllfi = 46.8 MeV. We see from Fig. that fx e ^ 
in the gCFL phase, indicating a nonzero electron den- 
sity. As discussed in detail in Refs. H, @, the negative 
Q-charge of the electrons is balanced by that of unpaired 
bu quarks, which have Q = +1, occurring in a narrow 
shell in momentum space. There is a larger shell in mo- 
mentum space, whose width grows with increasing M 2 / fi 
in the gCFL phase, within which there are unpaired fed- 
quarks. This "blocking region" of momentum space does 
not contribute in the Ai gap equation, and Ai is conse- 
quently driven down. The gap equations and neutrality 
conditions are all coupled, and the consequences of the 
reduction in Ai and the increase in fi e are manifest in all 
the curves in Figs. and |3] 

As Mg/n increases further, the gCFL solution even- 
tually ceases to exist at M 2 Jfi = 139 MeV. The gCFL 
solution to the gap equations is a minimum of Q with 
respect to variation of the A's for M^/fi < 139 MeV, be- 
comes an inflection point at M^/fi = 139 MeV, and for 
larger Ml/ fx there is no such solution. This is analyzed 
in greater detail in Refs. 0,0- The fact that the gCFL 
solution disappears indicates that there should be some 
other minimum with lower free energy, and indeed as 
shown in Fig. 0] we find that a first-order phase transition 
at which the gCFL free energy crosses above that of un- 



FIG. 3: Chemical potentials fi e , ^3 and /ig as a function of 
Mllfi at T = 0, for the (g)CFL and (g)2SC solutions of Fig.H 
(fi3 = in the (g)2SC solution and so is not shown.) Beyond 
the Ml / fx at which the g2SC solution ends, the dashed curve 
shows the chemical potentials (/i e 7^ and fis = 0) for neutral 
unpaired quark matter. 




50 100 150 

M^/jU [MeV] 



FIG. 4: Free energies of the (g)CFL and (g)2SC solutions of 
Fig. |21 at T = 0, relative to that of neutral unpaired quark 
matter. 

paired quark matter has occurred at M^/fi = 125 MeV, 
indicated in Fig.[21by the vertical line. In Ref. [||, the first 
order phase transition and the termination of the gCFL 
phase at a point of inflection of the free energy occur at 
M 2 s j [i = 130 MeV and M 2 J[i = 144 MeV respectively. 
Therefore, the errors in these quantities introduced by 
the small- Ml / /j, 2 approximation, used in Refs. but 
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not here, are about 4% at these values of M 8 2 //j. 

As in Ref. 0, we find an additional neutral 2SC so- 
lution, whose gap parameter A3 and free energy are 
shown in Figs. [3] and H At M s = 0, the 2SC gap is 
2 1 / 3 A [1113. For M s 2 //x below the 2SC^g2SC tran- 
sition at Mg/fj, = 112 MeV, rd-gu and ru-gd pairing 
occur at all mome nta* above this transition, in thegap- 
less 2SC phase [35|, tMl , there is a blocking region [33 in 
momentum space in which one finds unpaired rd and gd 
quarks, and A3 drops precipitously. An analogue of the 
g2SC pha se |3S| (and in fact an analogue of the gCFL 
phase 39]) were first analyzed in contexts in which they 
were metastable, but it was shown in Ref. [35| that the 
g2SC phase could be stabilized in two-flavor quark mat- 
ter by the constraints of neutrality. However, we see in 
Fig- IH that in this three- flavor quark matter setting, the 
(g)2SC solution everywhere has a larger free energy than 
the (g)CFL solution at the same Mj//x, and is therefore 
metastable. The value of M%//j, at which A3 — > in the 
g2SC solution is less than 1 MeV below the Mjf//x at 
which the gCFL phase becomes metastable. In contrast, 
in Ref. @ the g2SC solution persists to an Mf/fj, that 
is less than 1 MeV above that at which the gCFL phase 
free energy crosses that of unpaired quark matter. This is 
the one instance where the small- M s 2 //j 2 approximation 
made in Ref. 6] leads one (slightly) astray, as it pre- 
dicts a (very narrow) M|//x- window in which the g2SC 
phase is favored and we find no such window. However, 
we shall see below that the physics at values of M^/fj, 
that are this large compared to Ao is anyway not robust, 
changing qualitatively with increasing Ao. 

We now investigate how our zero temperature results 
change if we vary the strength of the coupling, and hence 
Ao. In Fig. EI we show the gap parameters as a function 
of Mj/// with A = 40 MeV. We have changed the scale 
on both the horizontal and vertical axes by the same 
factor of 40/25. We see that the CFL^gCFL transition 
again occurs at M^/fj, ~ 2Ai, and that the shape of the 
curves in the gCFL region is qualitatively as before, when 
suitably rescaled. However, at large values of M^/fi we 
now find a g2SC window: the gCFL free energy crosses 
above that of the g2SC phase — whose gap parameter A 3 
is also shown in Fig. [5] — at the vertical line in the figure, 
and the g2SC gap vanishes only at a larger Mf/fi. (We 
have located the vertical line by comparing free energies, 
as we did in Fig. 0] but we shall not give the figure.) 

If we reduce Ao from 25 MeV, rather than increasing 
it, and rescale both axes of Fig. [5] by the same factor 
by which we reduce Ao, we obtain a figure that looks 
qualitatively like Fig. [21 We conclude that stronger in- 
teraction tends to favor a g2SC window at large values 
of Mg/fi, whereas weaker interaction disfavors it. The 
boundary between the two cases is at Ao = 25 MeV in 
our model. 

It is interesting to ask what happens at still larger 
Ao. We show the gap parameters in our model with 
A = 100 MeV in Fig. El We see the CFL^gCFL tran- 
sition at Mg j /j, = 2Ai once again. The physics at and 
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M 2 s /jU [MeV] 

FIG. 5: Gap parameters for a stronger interaction than that 
in Fig. HI chosen such that A = 40 MeV for M s = 0. As in 
Fig.|2| the CFL— +gCFL transition occurs where M s 2 /^ ~ 2Ai. 
Here, however, the gCFL phase becomes metastable (at the 
thin vertical line) at a value of M% / fj, above which there is a 
g2SC solution. There is a first order gCFL^g2SC transition 
at the thin vertical line, followed at a larger /fj, by a second 
order transition at which the g2SC gap A3 vanishes. 



beyond the large- M^/fj, boundary of the gCFL regime 
is now qualitatively different. This regime corresponds 
either to very large values of M s , or else to such small 
values of fj, that the hadronic phase will likely have taken 
over, making the right half of this plot somewhat aca- 
demic. One reason it is of interest, however, is simply the 
fact that we can draw it: had we made a small M^/fi 2 
approximation as in Refs. 0, |(J 0, this regime would 
be inaccessible. The figure shows a sequence of phases as 
Ml) [i is increased above the gCFL phase: (i) the gCFL 
phase ends at a second order phase transition at which 
Ai — ► 0, above which we find a uSC window in which 
both A2 and A3 remain nonzero; (ii) next, A2 — > at 
a second order phase transition at which the uSC phase 
is succeeded by the 2SC phase; (iii) finally, there is a 
2SC— >g2SC transition. This sequence of phases agrees 
with that found in Ref. [TBI ] at large Ao, in a calculation 
done using a small- / [i 2 approximation pushed beyond 
its regime of validity. Our present calculation can be ex- 
tended (within its model context) to arbitrarily large M s . 
Indeed, as can be seen in Fig. {7\ we find a transition to 
two-flavor quark matter, with zero strange quark density, 
at M s = 547 MeV, corresponding to M 2 J \i = 598 MeV. 
Below this M s , we find the g2SC phase with unpaired 
strange quarks. Above this value of M s , we have two- 
flavor g2SC quark matter and further increase in M s has 
no effect on the physics. 

As Ao is increased from 40 MeV to 100 MeV, going 
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FIG. 6: Gap parameters for a still stronger interaction, 
with A = 100 MeV for M B = 0. As in Fig. H the 
CFL— >gCFL transition occurs where Ml/fi ~ 2Ai. At 
Ms /(i = 375 MeV, Ai vanishes at a second order gCFL^uSC 
transition. Then, at M^/fj, = 405 MeV, A2 vanishes at a 
second order uSC^2SC transition. At M^/fi = 449 MeV, 
there is a second order 2SC— >g2SC transition. And finally, 
for Ml /(i > 598 MeV, corresponding to M s > 547 MeV, 
there are no more strange quarks present in the system, as 
shown explicitly in Fig. |7| and further increase in M 3 changes 
nothing. 

from Fig.OOto Fig. El the first qualitative change to occur 
is that the gCFL phase ends at a second order transition 
at which Ai — ► 0, instead of ending at a first order tran- 
sition. Above this Ao, the phase diagram includes a uSC 
window separated from the (g)2SC phase by a first or- 
der phase transition. At a somewhat larger Ao, this first 
order phase transition becomes second order. At a still 
larger Ao, the interaction is strong enough to have g2SC 
pairing in the two flavor quark matter that our model 
describes for M 8 — > 00, and the physics is as in Fig. H3 



IV. THE PHASE DIAGRAM AT NONZERO 
TEMPERATURE 

We now explore the solutions to the gap equations and 
neutrality conditions at nonzero temperatures. As in the 
previous section, we begin with a coupling chosen so that 
Ao = 25 MeV. The phase diagram for this value of the 
coupling is given in Fig. Q to which the reader should 
refer in this section. We constructed Fig. Q by first mak- 
ing plots of the gap parameters and chemical potentials 
versus j \x at many values of T, and versus T at many 
values of M^/fi. In this Section, we present and discuss 
several of these "sections" of Fig. ^ enough to under- 



FIG. 7: Number densities of quarks with flavors u, d and s 
(in each case summed over the three colors) as a function of 
Ml I All parameters are as in Fig.© with A = 100 MeV. 
We see that the number densities are equal only in the CFL 
phase, and see that for Mj /fj, > 598 MeV there are no strange 
quarks present. 



2SC Solution 




FIG. 8: Gap parameters as a function of Ml / > at T = 2 MeV, 
with all other parameters as in Fig. [5] 



stand the many features of the phase diagram. We then 
show phase diagrams for Ao = 40 and 100 MeV. 

We start by turning on a small temperature, T = 
2 MeV, and seeing how the plots of gaps and chemi- 
cal potentials, shown in Figs. Eland El change from those 
at zero temperature, Figs. El and El We see many inter- 
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FIG. 9: Chemical potentials as a function of Mg/fi at T = 
2 MeV, with all other parameters as in Fig.EI At small M s 2 //x, 
He and (13 are close to 



esting changes already at this relatively small tempera- 
ture. The CFL^gCFL transition seen at zero tempera- 
ture in Figs. |21 and has become completely smooth at 
T = 2 MeV: there is no sharp difference between CFL 
and gCFL at nonzero T. Furthermore, note that the zero 
temperature transition from 2SC to g2SC is also washed 
out. This makes sense: at T = 2 MeV, it makes no phys- 
ical difference whether a certain fermionic quasiparticle 
is gapless or has a gap that is nonzero but smaller than 
2 MeV. So, although in Fig.Qwe have shown the values of 
M 2 j [i where quark quasiparticles become gapless within 
the CFL, uSC and 2SC phases, these dashed lines have 
physical significance only where they intersect T = 0. 

We see in Fig. [5] that Ai vanishes at a second order 
CFL-mSC transition at M|//x = 85.1 MeV. Another 
way to say this is that for M^/fx < 85.1 MeV, the critical 
temperature at which Ai vanishes upon heating must be 
greater than 2 MeV, whereas for Af s 2 /// > 85.1 MeV, this 
critical temperature is less than 2 MeV. As at T = 0, 
there is a first order phase transition, denoted in Fig. |S] 
by a thin vertical line, but here it is a first order phase 
transition between the uSC phase at M^//i < 125 MeV 
and the 2SC phase at M 2 J > > 125 MeV. The 2SC phase 
ends at a second order phase transition where A3 — * at 
Ml j [i — 129 MeV. This means that there is a regime of 
Ml) 'n at which there is 2SC pairing at T = 2 MeV, but 
no pairing at T = 0. We investigate this further below. 

We see in Fig. [5] that at small M^/fi, the chemi- 
cal potentials /z e and /Z3 are both close to Mg/4fi at 
T = 2 MeV. This is qualitatively different than their 
zero temperature behavior /i e = ^3 = 0. In comparison, 
/i e takes on the value in unpaired quark matter. 

Thus, already at a temperature of only 2 MeV there is no 




FIG. 10: Gap parameters as a function of T at M^/fji = 
15 MeV. 
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FIG. 11: Chemical potentials as a function of T at M^/fj, = 
15 MeV. 



sense in which /i e is small. We see that at larger M 2 / fj,, 
H e and /X3 diverge as at zero temperature, but they do 
so smoothly and they diverge from Mj/4/i, rather than 
from 0. Since we have found that the CFL phase has 
become a metal already at T = 2 MeV, it is natural to 
ask at what temperature the insulator-metal crossover 
occurs. We answer this question at M^/fi = 15 MeV in 
Figs. El and ^2 The latter figure shows a rapid insu- 
lator to metal crossover occurring between T = 1 MeV 
and T = 2 MeV, with fi e , ^ and /ig all changing. We 
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FIG. 12: Gap parameters versus temperature at Mg/fx = 
70 MeV. 

shall discuss this crossover at length in Section It 
can be understood analytically, and we shall see in Sec- 
tion[V]that the reason that /j, e takes on the value Mg/4fi 
is quite different from that in unpaired quark matter. 
Above the crossover, fj, e changes little as T increases fur- 
ther but both /is and /is decrease in magnitude. This 
occurs because at larger temperatures the gap parame- 
ters decrease, as seen in Fig. El and as the gap parame- 
ters vanish color neutrality occurs with fx 3 = fj, 8 = o M 
whereas electrical neutrality still requires a nonzero fi e . 

We see in Fig. ^| that the gap parameters change lit- 
tle at the low temperatures at which the CFL phase is 
undergoing its insulator to metal crossover. Although it 
is not really visible in the figure, we find that A2 > Ai 
for < T < 6.28 MeV, and A x > A 2 at higher tempera- 
tures. At T = 16.46 MeV, A2 vanishes at a second order 
phase transition and we find the dSC phase. Then, at 
T = 16.81 MeV, Ai vanishes, yielding the 2SC phase. 
The final phase transition, at which A3 vanishes, occurs 
at T — 17.73 MeV. This ordering of phase transitions 
is in qualitative agreement with that found in Ref. 0] 
using a Ginzburg-Landau approximation, and is in dis- 
agreement with the results of Ref. > m which no dSC 
regime was found. In order to gain confidence in the ac- 
curacy of our calculation and in the existence of the dSC 
phase, in Sect ion IVIl wc make a detailed and quantitative 
comparison between our results and those obtained via 
the Ginzburg-Landau approximation. 

At larger values of M^/fi, the ordering of phase tran- 
sitions as a function of increasing temperature changes. 
For example, if we consider M^/fj, = 70 MeV, in the 
gCFL phase at T = with Ai < A 2 < A3, we see in 
Fig. El that as the temperature is increased, Ai van- 
ishes first, then A 2 and then A3, meaning that the phase 
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FIG. 13: Chemical potentials versus temperature at Mf//x = 
70 MeV. 

which intervenes between CFL and 2SC is uSC, not dSC. 
This order of phase transitions is unsurprising, given that 
Ai < A 2 < A 3 at T = in the gCFL phase. All three 
transitions are second order transitions in mean field the- 
ory. Fig. El shows the chemical potentials as a function 
of increasing temperature at M^/fi = 70 MeV. Whereas 
the T = CFL phase undergoes an insulator to metal 
crossover as it is heated, the gCFL phase is already a 
metal at T = 0. 

In both Figs. El an d El we see a sequence of three 
second order phase transitions. The first of these, a 
transition from the CFL phase to either the dSC or the 
uSC phase, is likely not significantly affected by gauge 
field fluctuations, because the same gauge symmetries are 
unbroken (the U{1)q symmetry) and broken (the other 
eight gauge symmetries) on both sides of the transition. 
It is therefore an interesting question for future work to 
consider the order parameter fluctuations at this transi- 
tion, asking whether they render it first order or, if not, 
determining its universality class. The two mean field 
transitions occurring at higher temperatures in Figs. ITHI 
and El will be qualitatively affected by gauge field fluc- 
tuations, as at each of them there are gauge symmetries 
that are broken on the low temperature side of the tran- 
sition and restored above the transition. Gauge field fluc- 
tuations will presumably make these transitions first or- 
der, and shift their critical temperatures upward. These 
effects will be significant, because the relevant gauge 
fields are strongly coupled 0, H3 • 

We have found that at small M^/fj,, A 2 vanishes at 
a lower temperature than Ai whereas at larger / n, 
these two transitions occur in the opposite order. There 
must, therefore, be some Mj//x at which both vanish at 
the same temperature. We see in Fig.^that this "doubly 



14 



% 20 



03 



o 







- 






- 


_ 






_ 




|2SC 




A, A 2 



5 10 15 

T [MeV] 



FIG. 14: Gap parameters versus temperature at M^/fi = 
130 MeV. At this large value of M^/li, both Ai and A 2 are 
zero at all temperatures. A3 = at T = 0, but this gap 
parameter is nonzero for 2.37 MeV < T < 13.08 MeV. 



critical" point occurs at T — 15.3 MeV and M^/fi — 
29.4 MeV. 

Having followed what happens upon heating the CFL 
phase at M^/fx = 15 MeV and upon heating the gCFL 
phase at M^/fx = 70 MeV, in Fig. [21 we consider heat- 
ing quark matter at M% / fx = 130 MeV, which is unpaired 
at T = 0. We see that A3 becomes nonzero at a second 
order phase transition, and then vanishes at a higher tem- 
perature at a second second order phase transition. This 
behavior has been described previously |35l |40| , and can 
be understood as follows. At T — 0, the u and d Fermi 
surfaces in the unpaired quark matter are too far apart to 
allow 2SC pairing. However, as we increase T, we excite 
u quarks above the u Fermi surface, and d holes below the 
d Fermi surface. This smearing of the separated Fermi 
surfaces assists u-d pairing, and A3 turns on at a nonzero 
temperature. Of course, at a still higher temperature the 
A 3 condensate melts. 

The final slice of the phase diagram of Fig. ^ that we 
shall show explicitly is a plot of the gap parameters as 
a function of M^/ll at T = 4 MeV, shown in Fig. US] 
By comparing this figure with Fig. [HI we see qualita- 
tive changes in the physics between T = 2 MeV and 
T = 4 MeV: the 2SC phase now takes over from the 
uSC phase not via a first order phase transition, but in- 
stead via a second order phase transition at = 
123.2 MeV at which A2 — > 0. This means that a line 
of first order phase transitions, present at lower temper- 
atures, has turned into a second order transition at a 
tricritical point. This tricritical point is shown by a dia- 
mond in Fig.^ and is located between T = 3.9 MeV and 
T = 4.0 MeV. 




FIG. 15: Gap parameters versus M^/fx at T = 4 MeV. 



The results that we have presented up to this point in 
this section, plotted in Figs. 0D through ^[ constitute a 
description of all of the features at T 7^ depicted in 
the phase diagram given in Fig.^ We now ask how this 
phase diagram changes as we vary the strength of the 
interaction, and hence Ao. If we reduce Ao, there are no 
qualitative changes as long as we rescale the vertical and 
horizontal axes of the phase diagram by the same factor 
that we reduce Ao. As at T = 0, however, increasing 
A leads to qualitative changes in the phase diagram at 
large / fx, indicating that the details of the large Ml j ll 
regions of the phase diagram are not robust predictions 
of our model. 

With Ao = 40 MeV, in Fig. we see the same three 
special points as in Fig. ^ a quantum critical point sep- 
arating the CFL and gCFL phases at T = 0, a doubly 
critical point at which the Ai — > and A2 — > transi- 
tions cross, and a tricritical point at large /fx at which 
a line of first order phase transitions ends. Since both 
axes of Fig. 1161 have been rescaled by 40/25 relative to 
Fig. [2 the two fi gures are qualitatively similar: the one 
qualitative change occurs at large Ml /ll, where the g2SC 
phase extends down to T = 0, as we have already seen in 
Fig- El (At T = 0, there is a sharp distinction between 
2SC and g2SC, and in this instance the phase is g2SC.) 
The most interesting quantitative change is a change in 
the slopes of the Ai — > and A2 — > transitions on the 
phase diagram at small M|//x, which pushes the dou- 
bly critical point somewhat down in temperature. This 
effect is more clearly visible at stronger coupling, with 
A = 100 MeV, as shown in Fig. [T7| We shall explain 
this quantitatively in Section IVII 

With A = 100 MeV, in Fig. 1X71 we see further qual- 
itative changes at large M^/li, corresponding to those 
at T = shown in Fig. [SJ Now, the uSC phase and 
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FIG. 16: Phase diagram of dense neutral quark matter in the (Mf / fi, T) plane, with a coupling chosen such that Ao = 40 MeV. 
(All parameters except Ao are the same as in Fig. Only phase transitions are shown — the analogues of the dashed and 
dotted lines in Fig Q are not given. 



the 2SC phase both extend to T = 0. And, at T = 
there is a 2SC regime separated from g2SC by a quantum 
critical point, like that separating the CFL and gCFL 
phases at T = 0. If we start at T = in the g2SC 
phase and heat the system, A3 at first increases, before 
decreasing at higher temperature, eventually vanishing 
at the upper phase transition shown in Fig. 1171 If we 
extended the phase diagram to M^f/j, — > 00, this crit- 
ical temperature would become Ms-independent in the 
limit. The small remaining M s -dependence at the largest 
Mg/fx we show is easily understood: even though there 
are no strange quarks present at T — at these large val- 
ues of M s , meaning that the T = physics has become 
M s -independent, strange quarks can still be excited at 
nonzero temperature. 

As Ao is increased from 40 MeV to 100 MeV, the phase 
diagram changes continuously from that of Fig. llfil to that 
of Fig. El First, the uSC phase reaches the T = axis. 
Next, the tricritical point indicated by the diamond in 
Fig. ^retreats down to T = 0. All the while, the (g)2SC 
region is extending farther and farther to the right, even- 
tually to M s — > 00 when the coupling is strong enough to 
allow ud pairing even once there are no strange quarks 



present. 

We shall discuss the implications of the phase diagrams 
that we have found in Section VII. First, however, our 
investigation has raised several interesting questions that 
we have been able to address analytically, as we describe 
in Sections |V| and EH 



V. HEATING THE CFL PHASE: 
UNDERSTANDING THE INSULATOR TO 
METAL CROSSOVER 

We have seen in Fig. ^] that as CFL quark matter is 
heated, it undergoes a crossover from an insulator, with 
H e exponentially small, to a metal, with fi e ~ M^/Afi. 
With A = 25 MeV and M s 2 / 'fi = 15 MeV as in 
Fig. ^] the crossover occurs between T = 1 MeV and 
T = 2 MeV. This insulator to metal crossover has been 
seen previously in Ref. |l5j ; our goal here is to understand 
it analytically. 

We see from Fig. E3 that the gap parameters change 
little between T = 1 MeV and T = 2 MeV, whereas 
the chemical potentials change dramatically. We shall 
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FIG. 17: Phase diagram of dense neutral quark matter in the (Mf / fj,, T) plane, with a coupling chosen such that Ao = 100 MeV. 
(All parameters except Ao are the same as in Figs. ITIand llOI 'l The two triangles on the T = axis indicate quantum critical 
points at which modes become gapless. The first, which separates the CFL and gCFL phases, is familiar. The second separates 
2SC and g2SC phases, as shown in Fig. [S] 



explain the variation of the chemical potentials, treat- 
ing the gap parameters as T-independent. We note from 
Fig. II II that during the crossover, [x e increases from near 
to near while ^ increases by the same amount 

and us increases by half as much. This tells us that Q 
and (JHJ are maintained, and hence it is the combination 
Mq °f ® that is changing. We recall that at T = the 
contribution of the quark matter to the free energy Q, is 
independent of /j,q, and this "plateau" is only curved by 
the small contribution of the electrons to the free energy, 
of order which favors /i e = 0. Above the crossover we 
find /i e ~ which is on the plateau but away from 

the point on the plateau favored by the electron neu- 
trality condition. In order to understand the crossover, 
then, there must be thermally excited Q-charged quasi- 
particles whose neutrality condition favors fi e ^ 0, and 
we must see the curvature of the free energy plateau due 
to these quasiparticles "take over" from that due to the 
electrons. We hrst show that the CFL quark quasiparti- 
cle excitations have the desired effect, and then in Sec- 
tion IV Al we show that the thermally excited charged 
mesonic pseudo-Nambu-Goldstone boson excitations of 



the CFL phase play a negligible role. 

The quarks with nonzero Q are the rs and bu which 
pair with gap parameter A2, and the rd and gu which 
pair with gap parameter A3. This means that in the 
CFL phase there are two quasiquarks with Q = +1 (one 
a linear combination of bu quarks and rs holes; the other 
a linear combination of gu quarks and rd holes) and two 
quasiquarks with Q = — 1 (one a linear combination of 
rs quarks and bu holes; the other a linear combination of 
rd quarks and gu holes). 

We now evaluate the dispersion relations of these ex- 
citations, and estimate the number density of thermally 
excited charged quasiquarks. In this Section, we shall 
follow Refs. (Mill and include the nonzero strange quark 
mass only via its effect as a shift in the chemical poten- 
tials of the strange quarks. We have seen in Section ITTT1 
that this is a good approximation in the CFL phase. The 
dispersion relations of these four quasiparticles are given 
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by 



quasiparticles to the Q-charge density is then 



. 1 ( M 2 S 

£(rs-bu)[P) — ±2 I Mrs 7^ ~ Hbu 



£(rd-gu){p) = ±|(A*9« - Mrd) + J (p ~ ^(rd-gu)) 2 + A§ 

(22) 

where pbu-ra = (pbu + Mrs)/2 and p, r d-gu = (/Mi + AV)/ 2 - 
Upon substituting the definitions (|19l) and the relations 
Q and which are maintained through the crossover, 
these become 



2/i 



£(rs-bu)(p) = ± \L l e-^ L )+\l[p-H+-^ L ) + A2 



6/i 



ff(rrf- s „) (p) = ±Me + y - M + + A| . (23) 

At the temperatures of interest, these excitation energies 
are all greater than T, and so only the lowest energy exci- 
tation with each Q charge matters — the number density 
of the higher energy excitations is exponentially smaller. 
Labeling the quasiparticles by their Q charge, the exci- 
tation energies of the lowest lying charged quasiquarks 
are 




corresponding to the (rs-bu) quasiquark with Q = +1 
and the (rd-gu) quasiquark with Q = — 1. We now eval- 
uate the T-dependent contribution of these quasiquarks 
to the free energy fl given by JT3J. Because the quasi- 
quark energies are much larger than the temperatures of 
interest, the Boltzmann factors are small, the integral is 
dominated by p near the minimum of e(p), and we can 
use the saddle-point approximation. We find that, for 
example, the contribution of the Q = +1 quasiquark to 
fl is given by 



7T ./o 



dp p 7Tn(l + e 



-\e +1 (p)\/T\ 



-1 / dpp 2 Te^^\l T 
71-2 Jo 



(25) 



where /i = p — /6/i. The contribution of the Q = — 1 
quasiparticle is analogous. The contribution of these two 



~ = \^ArM 2 e - ( A +Me -^/2 M )/T ) (26) 

(Q = -1) - ^^! e -(A-,e)/T (2?) 



We have set A3 = A2 = A, a good approximation 
throughout the crossover as shown in Fig. ^| Q- 
neutrality is a balance between the charge densities of 
these two quark quasiparticles and the electrons. 

We can now see that what drives the system to fi e 7^ 
is the fact that the lightest Q = +1 and Q = — 1 
quark quasiparticles have dispersion relations with dif- 
ferent gaps when M s ^ 0, and consequently contribute 
charge densities of different magnitudes when thermally 
excited. If we attempt to set fi e — at T ^ 0, there 
are more Q = +1 quasiparticles present than Q = — 1 
quasiparticles, and the system is not neutral. To achieve 
neutrality, fi e must be increased as this increases the den- 
sity of Q = —1 quasiparticles, decreases the density of 
Q = +1 quasiparticles, and adds electrons, which have 
Q = — 1. This is described by the neutrality condition 



2V2irATfl 2 



PeT 2 



which we can solve for pi e (T) if we take A to be T- 
indepcndcnt. 

Let us now investigate the implications of this result. 
At very small T, the quark quasiparticles are exponen- 
tially rare, and those with Q = — 1 are exponentially 
rarer than those with Q = +1. The Q = +1 quasipar- 
ticle density is balanced by the electron density, and /i e 
is exponentially small. However, the quasiparticle densi- 
ties are proportional to /1 2 whereas the electron density 
is not. This means that at the T at which /i e starts 
to take off, the quasiparticle Boltzmann factors are still 
rather small. Once T is large enough that /i e approaches 
M^/4p, however, even though the individual Boltzmann 
factors continue to rise rapidly as T increases further, 
the sinh factor in l|28l) becomes small. Neutrality at this 
point is primarily a balance between the densities of the 
Q = +1 and Q = — 1 quasiparticles, with electrons can- 
celling only the small difference between their densities. 
The result, seen already in Fig. Illland shown in greater 
detail in Fig. 1181 is a crossover in which fi e is at first 
exponentially small, then rises rapidly, and then satu- 
rates as it approaches M 2 /4/j,. We see in Fig. 1181 that 
the equation 1)281 for [i e that we have derived , making 
approximations as described, gives a very good descrip- 
tion of the numerical solution of the full coupled gap and 
neutrality equations. This demonstrates that (|28l) pro- 
vides us with a good analytic description of the insulator 
to metal crossover that CFL quark matter experiences 
when heated. 
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FIG. 18: Comparison between the analytic results (dashed 
curves) and the numerical results (solid curves) for as a 
function of temperature with Ao = 25 MeV for several values 
of M a in the CFL phase. The analytic results were obtained 
from 12811 , and the numerical results were obtained by solving 
the full coupled gap equations and neutrality conditions, as in 
previous Sections. In evaluating /i e from 128H we have taken 
A to be the average of A2 and A3 at the midpoint of the 
crossover. 



We have set the electron mass to zero in Fig. ^| and 
throughout. Including it means that it takes a larger /i e 
to achieve a given electron density, pushing all the curves 
in Fig. 1181 very slightly upwards. With an electron mass 
as in nature, the effect on the curves is invisible on the 
scale of the plot. 



A. Contribution of charged mesons 

As described in Section Til Bl there are other charged 
excitations in the CFL phase. Among the pseudo- 
Nambu-Goldstone bosons, there are mesons with the 
quantum numbers of the ^ and K . We have neglected 
the contribution of thermally excited charged mesons to 
the charge density in the derivation of 128|) . We now 
investigate this approximation. 

The dispersion relations of the charged mesons, to- 
gether with that for the i^°-meson which we shall also 




T [MeV] 

FIG. 19: Comparison between the analytic estimate for 
versus T with and without charged meson contributions. The 
solid curve, with no meson contributions, is the same as the 
M s = 100 MeV analytic estimate in Fig. EU The dashed 
curves include the contribution of thermally excited charged 
pions and kaons, for three values of the parameter \ described 
in the text that parameterizes the instanton contribution to 
the meson masses. 

need below, are given by [H |H H HI |H |53 

e^± (p) = ±fi e + ^Jv 2 p 2 +M 2 ± 

M 2 I 

e K ± (p) =iftT^ + ^v 2 p 2 + M 2 K± 

M 2 1 

e K o(p) = -^ + ^v 2 p 2 + M 2 KO , (29) 

where v 2 = 1/3 at high density HU . The meson masses 
in the CFL phase are given by 

M 2 ± = a(M u + M d )M s + X (M U + M d ) 
M|± = a(M u + M s )M d + X {M U + M s ) 
M 2 R0 = a(M d + M S )M U + X (M d + M s ) . (30) 

Here, a = 3A 2 /ir 2 f 2 with f 2 = (21 - 8 log 2)^ 2 /36tt 2 
at high density [2f|, which yields a = 0.0175 for A = 
25 MeV at n = 500 MeV. And, \ parameterizes the 
contribution of [/(l^-breaking instanton effects which 
generate (qq) condensates and therefore contributions to 
meson masses in the CFL phase H?}, via the 
't Hooft interaction which contributes couplings of the 
form A 2 (qq). The magnitude of x is not well known, as it 
depends on the instanton size distribution and instanton 
form factors at nonzero density. It has been estimated 
to lie in the range 1 MeV < x < 100 MeV p l l27 | . In 
Fig. ^3 we include the contribution of all four charged 
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mesons to the Q-charge neutrality condition, determin- 
ing the density of thermally excited bosonic quasiparti- 
cles from the dispersion relations. The contribution to 
the charge density from the mesons is 

2^1 dp p2 ( exp [e K+ (p)/T] - 1 ~ cxp [e K - (p)/T] - 1 ) ' 

(31) 

and that from the pions is analogous. We then solve for 
He vs. T, taking M s = 100 MeV and M u = 5 MeV and 
M d = 10 MeV. We plot the results for \ = 1, 2 and 5 
MeV in Fig. 1151 

Adding the charged mesons adds new charge carriers, 
and so the simplest expectation for their effects is that 
a smaller ji e will be required in order to neutralize the 
imbalance in the fermionic quasiparticle sector. We see 
in Fig. El that this expectation is borne out above the 
crossover, but not below. Above the crossover, we see 
that n e is reduced relative to the results we obtained in 
Fig. El where we neglected the mesons. The mesonic 
contributions get less significant at larger values of x, 
as the mesons get heavier. They are already small for 
X = 1 MeV and are negligible by x = 5 MeV. At low 
temperatures, and for the smallest values of \, there is 
another effect to be understood. The dispersion relations 
for the K ± in (|29[1 indicate that the K + is easier to ex- 
cite than the K~ . This means that they behave like the 
quasiquarks, in the sense that at a nonzero temperature 
they contribute a positive Q-charge density, which must 
be cancelled. If x is very small, the K + charge density 
becomes significant at such a low temperature that its 
contribution can only be cancelled by electrons — the 
Boltzmann factors for all other excitations are still pro- 
hibitive. This means that if x is very small, [i e initially 
rises with temperature significantly more rapidly than in 
the absence of the mesons. We see this effect clearly in 
Fig. El f° r X — 1 MeV and still to a small degree for 
X = 2 MeV. For x = 5 MeV, fi e (T) is indistinguish- 
able at low temperatures from that in the absence of the 
mesons. 

We can summarize the results in Fig. El as follows. 
F° r X ~ 2 MeV, the thermally excited charged mesons 
have no significant effects at any temperature. But, if 
X ~ 1 MeV, which is at the bottom end of the estimated 
allowed range 1 MeV < x < 100 MeV |HE3, the K+ 
excitations contribute significantly at very low tempera- 
tures. 

At much larger temperatures, well above the insula- 
tor to metal crossover that we have analyzed, we come 
to the various critical temperatures that we have seen 
in Section Hvl and will analyze further in Section IVT1 At 
these temperatures, at which gap parameters vanish, it is 
well understood that the fermionic quasiparticles are the 
most important degrees of freedom. What we have shown 
is that at low temperatures also, the charged mesons are 
less important than the charged fermions among the ther- 
mal excitations, as long as x is not very small. 

Mesons can nevertheless play an important role if they 
condense [13 . At T = 0, the CFL phase is stable against 



meson condensation as long as Ek« > at p — 0. (K°- 
condensation yields the most stringent constraint.) This 
requires M|- D > corresponding to x ^ M%/4[j, 2 . 

For M 8 = 100 MeV, as in Fig. El this requires x > 
1 MeV. We are therefore justified in our neglect of K°- 
condensation in Fig. El If X > 3.6 MeV, there is no 
A"°-condensation in the CFL phase with Ao = 25 MeV 
for M s 2 / ' fji < 46.8 MeV, meaning that there is no K°- 
condensation at T = for all values of M 2 j /j, below the 
CFL— s-gCFL transition. (A°-condensation in the gCFL 
phase has yet to be analyzed.) It seems likely that x > 
3.6 MeV at accessible densities, since y is larger at lower 
densities and the analysis in Ref. [2J| which yields the 
estimate 1 MeV < x < 100 MeV becomes more reliable 
at higher densities. If the coupling is stronger, however, 
A"°-condensation becomes more likely. For example, for 
Ao = 100 MeV, there is no A°-condensation in the CFL 
phase only if x > 26 MeV. If A°-condensation were to 
occur, it delays the CFL^gCFL transition, increasing 
the Mg at which it occurs by a factor 4/3 if x = [l3[ . 
Further work remains to be done, for example extending 
the analysis of Ref. [R| to nonzero x an d to the gCFL 
phase. 

Another issue that remains to be investigated is the 
possibility of 7r~ condensation at nonzero temperature. 
We have seen that above the insulator to metal transi- 
tion, yit e w M s 2 /4/it. At zero temperature, this would lead 
to 7T~ condensation if x < Af s 4 /(16^ 2 (M U +Af d )), but this 
li e only arises at T / 0, and nonzero T acts to stabilize 
against meson condensation. 

Finally, all these issues should be investigated in an 
expanded model in which the 't Hooft interaction is in- 
cluded from the beginning in the free energy and hence in 
the gap and neutrality conditions, and the quark masses 
are also solved for dynamically. 



VI. THE GINZBURG-LANDAU 
APPROXIMATION 

In the M s 2 / fi — > limit, the three critical tempera- 
tures at which the gap parameters vanish become one. 
Near this critical temperature, where all gap parameters 
are small, a Ginzburg-Landau approximation can be em- 
ployed. The Ginzburg-Landau free energy was analyzed 
at M s = in Refs. |3, EJ, and was extended to small 
but nonzero M 2 / ' \i in Ref. 16]. We wish to compare 
our results at small but nonzero M 2 /fi to those obtained 
in the Ginzburg-Landau approximation, and to compare 
the coefficients in the Ginzburg-Landau potential (actu- 
ally, ratios of coefficients) in our model to those calcu- 
lated in QCD at asymptotic densities and thus at weak 
coupling in Refs. |3lll6|. 

The Ginzburg-Landau potential is parameterized 
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FIG. 20: Diagrams contributing to the Af and AiA| terms 
in the Ginzburg Landau potential 1321 , 



as [3 

n =a{A\ 
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(32) 



/? 2 (Af 



A^ 



where a = cvo(T — T c )/T c and where the coefficients e 
and ?7 are proportional to M 2 and /i e respectively and 
are therefore present only if M s ^ . The form of the e 
and r\ terms was derived in Ref. ,10] for QCD at asymp- 
totic densities. In the Ginzburg-Landau limit, in which 
all gap parameters are small, color neutrality occurs with 
Us and /is vanishingly small ^(|, and electrical neutral- 
ity requires a nonzero /Lt e , of order M 2 /4[i as in unpaired 
quark matter. At low temperatures, the most important 
consequences of M s and fi e (and also /13 and /ig) derive 
from the stress they put on pairing, as they (seek to) 
push Fermi momenta apart. Near T c , T ^> A. And, we 
are working at M 2 /fi <C T c , where the thermal smear- 
ing of the Fermi distributions is much greater than the 
splittings between Fermi momenta. In this regime, the 
effects of M 2 or /x e on the pairing of two quarks do not 
arise from the splitting between their Fermi momenta, as 
this is negligible. The effects of or fi e arise from the 
change in the average Fermi momenta of the two quarks, 
and hence in the density of states, that M 2 or /i e induces. 
For example, M 2 depresses the average Fermi momenta 
of u-s and d-s pairs, but does not affect u-d pairs. This 
explains the form of the e term in (|32[1 . On the other 
hand, fj, e > increases the average Fermi momenta of d- 
s pairs, while decreasing that of u-d and u-s pairs. This 
explains the form of the i] term. 

In Ref. Q it is shown that, in the weak coupling regime, 
Pi = fa- This result is valid beyond weak coupling, how- 
ever, as all that is required to demonstrate it is the pair- 
ing ansatz © and the mean field approximation. The 
quarks in the 2x2 blocks contribute through diagrams 
like the first in Fig. 1201 leading to a contribution propor- 
tional to Af + A| + A|. The quarks in the 3 x 3 block 
contribute through diagrams like both the first and sec- 
ond in Fig. 1201 leading to a contribution proportional to 
Af + A| + Af + 2A 2 A 2 + 2A 2 2 A 2 3 + 2A 2 A§. Adding all 
the diagrams, we conclude that fa = fa- 
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FIG. 21: Gap parameters squared versus T at M B = 30 MeV, 
namely M% /n = 1.8 MeV. This figure should be compared to 
the "schematic illustration" given in Fig. 2 of Ref. |lfij|. We 
see the three phase transitions separating the CFL, dSC, 2SC 
and unpaired phases. 



We now wish to confirm that (|32|l correctly describes 
the M s -dependent physics in our model. We do so by 
extracting the ratio 



c 



e 
V 



(33) 



from our results in three independent ways. If l|32|l is 
correct — that is if there are no M s -dependent terms 
missed — the three extractions of C should agree. We see 
in the phase diagrams of Figs.lTl[RTland[T7lthat the slopes 
of the three transition temperatures (i.e. dT c /d(M 2 /[i) 
at Mg/fj, = 0) are different, with that for the A3 — > 
critical temperature the shallowest and that for the 
A2 — * critical temperature the steepest. The ratios of 
these slopes can be extracted from the Ginzburg-Landau 
free energy £1 of 132|) and are given by 



1 : (6C - 5) : (6C + 4) 



(34) 



where we have used fa — fa- We can also read the slopes 
directly from our phase diagrams. The Ao = 25 MeV 
phase diagram of Fig. ^ yields the ratios 



1 : 10.1 : 19.1 



(35) 



which can be used to obtain two independent extractions 
of C, one from 1 : (6£ — 5) = 1 : 10.1 and the other from 
1 : (6£ + 4) = 1 : 19.1. The two are in perfect agreement, 
with both yielding ( = 2.52 for A = 25 MeV. 

A third extraction can be obtained from Fig. 1211 upon 
realizing that, according to l|32|) . in the CFL phase where 
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all three gap parameters are nonzero C is given by the 
ratio 



This ratio can be extracted at M 2 //j, — 1.8 MeV from 
Fig. |^ We have done this extraction at a number values 
of Mg I /i and fitted the results as 

^§ r ^ + 3 ^) 2 +1 , 2xl0 ^ff (37, 

which leads us to conclude for the third time that C = 
2.52 for A = 25 MeV. In ^3 we nav e extracted 
the M s 2 //i 2 and Mf/fi 4 corrections to the ratio in 
addition to extracting £. The coefficients we have ob- 
tained confirm that these higher order terms are negligi- 
ble as long as M 2 / fi 2 <C T c j ' \i or, equivalently given QJ, 
Mg I [i -C Ao . This indicates yet again that the Ginzburg- 
Landau free energy ((32(1 provides a good approximation 
to our results in the regime where it should be valid. 

Upon comparing our Fig. 1211 with the schematic illus- 
tration given in Fig. 2 of Ref . [0| , there can already be 
little doubt that the Ginzburg-Landau potential 1 (3*2) 1 cor- 
rectly describes the results that we have obtained by solv- 
ing the full gap and neutrality equations numerically, for 
T ~ T c and for M 2 / fi 2 < T c / fi. The agreement between 
our three extractions of £ makes this point quantitatively. 

We can also use Fig.[!""""]to check that our numerical re- 
sults are consistent with fa = fa. The Ginzburg-Landau 



potential (|32|l can first be used to show that in the CFL 
(and dSC) phases in Fig. |2] the slopes of the lines for the 
three (two) nonzero gap parameters are the same. Our 
results clearly satisfy this. Next, the Ginzburg-Landau 
potential (|32|) can be used to show that the ratio of the 
slopes of the gap parameters in the CFL phase in Fig. 1211 
to those in the dSC phase is (2 fa + fa)/ (3 fa + fa), and 
the ratio of those in the dSC phase to that in the 2SC 
phase is (fa + /?2 )/(2/3i + fa). These ratios cannot be 
extracted very accurately from Fig. 1211 given the nar- 
row windows within which the dSC and 2SC phases are 
found. However, they are consistent with 3/4 and 2/3, 
corresponding to fa — fa. 

In Ref. ^(|, the ratio £ is calculated using weak cou- 
pling methods, valid at asymptotic densities. The weak 
coupling result is £ = 2. From our numerical results, 
we have found Q = 2.52 at A = 25 MeV. We also 
find ( = 2.69 at A = 40 MeV and < = 3.87 at 
Ao = 100 MeV, extracting £ from the ratio (|36ll as in 
1)37(1. At weak coupling, C — 2 and the A3 and A 2 , lines 
are equidistant from the A 2 line in the CFL phase region 
of Fig. [21 As the coupling gets stronger, the A\ and A\ 
lines move downward/leftward, further away from the A 2 
line, and the ratio £ increases. 

Now that we are convinced that the e and 77 terms fully 
describe the M s -dependent physics in our model at small 
M s , we calculate e and rj, from diagrams like those in 
Fig. 03 After some calculation, the result is 
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2T c (p + n) 



cosh 
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(38) 



(39) 



In the weak-coupling limit, T c <C /i and the integrals are 
dominated by p within T c of fi. In this limit, it is easy to 
check that £ = e/77 = M 2 /2fifi e , which yields £ = 2 since 
fi e = Mg/4/i. At non-infinitesimal coupling, for example 
taking Ao = 25 MeV and reading the corresponding T c 
from Fig.Ql we can evaluate e and rj. We find £ = 2.55 if 
we use fj, e = M 2 /4/x, and if instead we obtain fj, e from our 
numerical results, we find £ = 2.52. Taking T c and /x e 
from our numerical results and evaluating £ using 



we find C = 2.52, 2.65, 3.84 for A = 25, 40, 100 MeV 
respectively. From our numerical results, we have found 
C = 2.52, 2.69, 3.87 at these values of Ao, extracting £ 
from the ratio (|36|l as in ((37(1 . The agreement between 
these determinations is a confirmation of the accuracy of 
our numerical methods. 

It is nice to see how quantitatively well the Ginzburg- 
Landau approximation describes the physics near T c in 
our model, as we have demonstrated. Furthermore, the 
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FIG. 22: Diagrams contributing to the e and rj terms in the 
Ginzburg Landau potential 1321 . 



value of one of the two ratios of coefficients that we have 
investigated, /9i//?2 = 1, is the same in our model and in 
QCD at asymptotic densities. And, the value of the other 
ratio C = e/r] is comparable in our model to its value in 
QCD at asymptotic densities for A = 25 and 40 MeV, 
becoming significantly larger only for quite strong cou- 
pling, as at A = 100 MeV. 

VII. IMPLICATIONS AND OPEN QUESTIONS 

The phase diagrams shown in Figs. HI HD and 1771 con- 
stitute the central results of this paper. They could be 
used in one of two different ways. If future theoreti- 
cal advances constrain the /^-dependent values of Ao and 
M s more tightly than at present, these phase diagrams 
(and those that interpolate between them for other val- 
ues of Ao) could be used to construct the phase diagram 
of nature. Or, if future astrophysical observations teach 
us that the phase diagram of nature must have certain 
features, for example must or must not include a cer- 
tain phase, then the phase diagrams we have constructed 
could be used to draw inferences about the magnitudes 
of A and M s . 

In thinking about the future phenomenological use of 
the phase diagrams that we have found, their complexity 
raises concerns. However, in most astrophysical contexts 
compact stars have temperatures much less than 1 MeV. 
At such low temperatures, which can be thought of as 
T = in a QCD context, the phase diagrams are more 
manageable. We have the CFL phase at asymptotic den- 
sities, with the gCFL phase taking over at lower densities, 
when /j, < M^/2A. If A lies at the large end of its es- 
timated range 10 MeV < A < 100 MeV, it seems likely 
that hadronic matter will take over from gCFL (or even 
from CFL) meaning that the complexities that we have 
found in our phase diagrams at larger Mj//i will likely 
be superseded by the transition to the hadronic phase. 

If Ao lies at the lower end of its allowed range, then 
in Fig. ^ we find a straightforward transition to "un- 
paired" quark matter as the density is decreased further 
and M s 2 / /j, increases beyond the gCFL window. If we ex- 
tend our pairing ansatz, we will certainly find some pair- 



ing in this regime at sufficiently low temperature. For 
example, perhaps weak pairing between quarks with the 
same flavor plays a role 0, , or pe rhap s it is the crys- 
talline color superconducting phase |37l |42| that takes 
over from gapless CFL at lower densities. Recent devel- 
opments |42j | make the crystalline color superconducting 
phase look like the most viable contender |y. 

The one astrophysical context in which the full com- 
plexity of the phase diagrams that we have analyzed must 
be faced head-on is the physics of the proto-neutron star 
formed during a supernova, and in particular of the prop- 
agation of neutrinos therein. Phenomena encoded in the 
phase diagrams that we have analyzed could ultimately 
result in observable consequences in the time-of-arrival 
distribution of the neutrinos detected from a future su- 
pernova [53, 12^. |43| . The phenomenological implications 
of the complexity of the phase diagrams (with many 
phase transition lines, the doubly critical point, the tri- 
critical point and the insulator to metal transition as the 
CFL phase is heated) will have to be thought through 
in this context. The analytic treatment of the insulator 
to metal crossover and the Ginzburg-Landau analysis of 
the physics near T c for small Mj / 1 [iT c that we have pre- 
sented in Sections and IVII could prove valuable in this 
context. 

Our analysis leaves many open avenues of investiga- 
tion that must still be followed to their conclusions. The 
effects if gauge field fluctuations must be investigated, as 
must those of K °-condensation in the gapless CFL phase. 
The effects of the 't Hooft interaction should be included, 
and the quark masses should be treated as dynamical 
condensates to be solved for, rather than as parameters. 
The pairing ansatz should be generalized, for example to 
allow for a comparison between the free energies of the 
gapless CFL and crystalline phases in three-flavor QCD. 

Finally, the stability of the gapless CFL phase needs 
to be investigated further, along the lines of Ref. [44| . 
These authors find that the g2SC phase is unstable, as 
is the 2SC phase near the 2SC— >g2SC transition. This 
instability could reflect the known instability of the g2SC 
phase with respect to a mixed phase of charged compo- 
nents 45], it could reflect an instability with respect to 
the crystalline color superconducting phase, or it could 
reflect the existence of some inhomogeneous phase yet 
to be discovered. We know that the gCFL phase is sta- 
ble with respect to mixed phases, except perhaps at the 
largest values of M% / // where it is found in our phase dia- 
grams ja, |4(j . It seems likely that at these large values of 
M s 2 / \i the gCFL phase may anyway be superseded by the 
crystalline phase This could be clarified either by 

a direct free energy comparison of these two phases, or 
by a stability analysis of the gCFL phase as in Ref. [44| • 
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